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ABSTRACT 

Starting from a field theoretical description of multicomponent anyons with mutual 
statistical interactions in the lowest Landau level, we construct a model of interacting 
chiral fields on the circle, with the energy spectrum characterized by a symmetric matrix 
g a p with nonnegative entries. Being represented in a free form, the model provides a field 
theoretical realization of (ideal) fractional exclusion statistics for particles with linear 
dispersion, with g a p as a statistics matrix. We derive the bosonized form of the model 
and discuss the relation to the effective low-energy description of the edge excitations for 
abelian fractional quantum Hall states in multilayer systems. 
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1. Introduction 

Recently a field theory of anyons confined to the lowest Landau level (LLL) of a strong 
magnetic field was formulated [[|] . The authors of Ref . p] next mapped their field theory 
(which is effectively one dimensional) onto a circle obtaining a theory of a chiral fermion 
field with linear dispersion. The problem involves a (statistics) parameter g entering the 
wave function of anyons in the LLL ^ = Ili<i( z j — z j) 9 ^+{{ z i\) e ~ 1 '^~ Ji where $ + ({zj}) 
is a symmetric function of its arguments. The requirement that the function $ + should 
be a polynomial for all the states in the LLL implies that g G [0, 2). Remarkably, the field 
theory on the circle remains consistent in a wider interval of range of the parameter g, 
for all g > 0. The motivation to consider such an extension of the theory is the similarity 
of the above wave function with g odd to Laughlin wave functions of fractional quantum 
Hall (FQH) states 0. It suggests that for g odd, the model of a chiral field discussed in 
|I[ may be used to describe edge excitations for Laughlin states. Indeed, a bosonization 
of the theory enables one to recover correlation exponents as well as the low temperature 
thermodynamics |lj coinciding with those obtained with the use of the (chiral) Luttinger 
liquid description of edge excitations |5| (for a recent review, see M). 

In this paper we generalize the model in to the case of several species of particles. In 
Section 2, we discuss a field theoretical description of anyons of several species in the LLL 
described by the generic statistics matrix g a/3 (with off-diagonal elements corresponding to 
mutual statistical interactions ||), and then map the system onto a circle. Upon extending 
the possible values of the statistics parameters, we arrive at a model of interacting chiral 
fields, with the energy spectrum characterized by a symmetric matrix g a p with non- 
negative entries. As in the case of a single species, the theory may be reformulated in a 
free form (in terms of field operators obeying more complicated commutation relations). 
We discuss the construction of the corresponding Fock space. In Section 3 we show that the 
structure which is behind the above free field formulation is an ideal fractional exclusion 
statistics || [7|, |], ||, [10], [11], in the thermodynamic limit the model is equivalent to 
a system of free chiral particles with linear dispersion, characterized by the (bosonic) 
statistics matrix g a/3 . 
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It was observed in Ref. [EJ that as a result of the dimensional reduction, the system 
of two anyons in the LLL admits a description in terms of fractional statistics in one 
dimension. The latter statistics was originally introduced by Leinaas and Myrheim in 
the Heisenberg quantization of identical particles |TJ]]. It was also argued in Ref. |TJ| 
that the system of noninteracting particles obeying fractional statistics in one dimension 
is equivalent to a system of particles with long range inverse square interaction (see also 
flT5|). Using the asymptotic Bethe ansatz for the latter model [IB], an equation for the 
single state distribution function for fractional statistics in one dimension was derived 
(note that in spite of being derived from the one- dimensional model, the resulting single 
state statistical distribution does not make any reference to the space dimension) JF]. The 
same distribution function was also derived [§, || from the generalized exclusion principle 
and the statistics is now referred to as fractional exclusion statistics. The system of 
anyons in the LLL is a realization of this statistics, as can be seen from the equation of 
state and thermodynamic quantities || [K| || . 

The above explains the appearance of exclusion statistics in the model of chiral fields on 
the circle which has the same quantum numbers as the system of anyons in the LLL. Using 
the exclusion statistics representation of the model, we calculate the low temperature heat 
capacity. The connection with exclusion statistics enables one to describe the model in 
terms of equations having the form of asymptotic Bethe ansatz equations |16j. This 
implies a simple dynamics encoded in two-body scattering phases of a step-wise form. 

We discuss some peculiarities in calculating quantities like the dressed charge matrix 
due to a step-wise character of the two-body scattering phases. 

In Section 4 we consider the relation to the effective low-energy description of edge 



excitations in abelian multilayer FQH states |17], ||. In the bosonized form, our model 
describes a set of chiral boson fields all propagating in the same direction with the same 
velocity. For the case of odd diagonal entries and integer off-diagonal entries, the statistics 
matrix g a p can then be identified with the topological matrix K a p describing the edge 
excitations. This in particular implies an interpretation of edge excitations in terms 
of exclusion statistics and is consistent with the results of Fukui and Kawakami who 
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achieved the same identification between g a p and K a p for edge excitations in hierarchical 



FQH states 12 



2. Multicomponent field theory on the circle 

Consider a system of M species of anyons in the LLL, described by the generic statistics 
matrix g a p (a, ft = 1, . . . M), with N a the number of particles of species a. The wave 
function of the system has the form || 

* = n n^r - z i) 9aa n nc? - (A e, wi a . (i) 



Here zf = \JeB/2(x c * + iyf) is a dimensionless complex coordinate of the i th particle 
of species a (the charges and masses of all the species of particles are assumed to be 
the same). The symmetric matrix g a p is responsible for the statistics of anyons, with 
non-diagonal elements corresponding to mutual statistical interactions between particles 
of distinct species. The function <3> + is single-valued and symmetric with respect to an 
interchange of coordinates of the same species, arbitrary when coordinates of distinct 
species are interchanged (bosonic representation). 

The wave function ([!]) acquires a phase exp(i7rg aQ ,) under an (anticlockwise) inter- 
change of two particles of species a and a phase ex.Y>{i2itg a p) under winding a particle of 
species a around a particle of species ft for a ^ ft (provided that each of the above closed 
paths encloses no other particles). In order for the function $ + to be a polynomial for 
all the states in the LLL, one should choose the intervals g aa G [0, 2) and g a p G [0, 1) for 
a 7^ ft. g a p = and g a p = 5 a p correspond to sets of bosons and fermions, respectively. 

It will be convenient for us to work in the fermionic representation. The latter is 
obtained by the change g a p = 5 a /3 + X a /3 and absorbing a Slater determinant corresponding 
to the Sap-p&rt of g a p in (HD in the definition of the function $. This converts $ + into 
the function $_({zf }, {zf }, ...) which is antisymmetric with respect to an interchange of 
coordinates within the same species of particles. 

The angular momentum for the states ([ID has eigenvalues 



E 



i=i 2 /¥<* 



(2) 



where v? a are nonnegative integers, distinct within each species, corresponding to the 
fermionic part of the angular momentum associated with the function 

Following Ref. 0, we introduce the fermionic field operators <p a (z) obeying the fol- 
lowing anticommutation relations, 

{ ( Pa{z),<Pfi(*f)} = 0, &a(z'),<pl(z)} = 5 aP -e z ' z . (3) 
On expanding <p a (z) in powers of the angular momentum eigenstates 

OO 1 

iPa(z)='£- T =ipf^, (4) 
1=0 V 7Ti! 

these imply the usual fermionic anticommutation relations for the components: 

{<pf, rf} = {^, ^} = o , K, rft} = . (5) 
In terms of the field operators @, the angular momentum can be written as 



z9 * + \ H A "/3^ 
2 /3 



yj a (z) , (6) 



a J 

with the representation of the particle number operators as 

N a = J d 2 ze-- zz vi(z) Va (z) (7) 

(d 2 z = dxdy). By using the expansion (J3|) and the explicit construction of the fermionic 
Fock space generated by the operators fl5|), one can verify that the formula (D) recovers 
the correct angular momentum eigenvalues (§) for states with fixed particle numbers. 
Upon adding a harmonic potential (of frequency to), the Hamiltonian becomes 

H = ^cufJ2N a + aL, (8) 



where oof 1 = Ju>% + 4u 2 , u c = eB/m is the cyclotron frequency, and a = uf 1 — u c . To 
leading order in u 2 /u) 2 , which is assumed hereafter, 



u;^ ff ~u; c , a~u)' 2 /uj c , lu<^lu c . (9) 
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The wave functions ([!]) for g aa odd integers and g a p integers have the form of the wave 
functions describing FQH abelian states in multilayer systems, with zf playing the role 
of the complex coordinate of the 2 th electron in the a th layer (such wave functions were 
first discussed by Halperin |I8|). Motivated by this similarity, we consider an analytic 



continuation of the solutions for anyons in the LLL, allowing the parameters g a p to take 
any nonnegative values. For a single layer, the edge excitations are generated when the 
function $ + in (|I|) belongs to the space of symmetric polynomials in {zi} [Q, |J, which is 
relevant to the polynomial character of the functions $ + in (Q). 

We now consider a mapping of the system onto a circle. Using the operators (pf from 
(H), we define the fields on the circle (parametrized with an angular coordinate 9) by 

00 1 

xM = E -7trrtJ nae - (io) 

The commutation relations (|5]) then imply 

{x*(0),Xp(0')} = {xU0),4(6')} = 0, 

1 oo 
Z7F n=0 

where S^ eT (9 — 6') is the positive frequency part of the periodic 5-function. 

We introduce the Hamiltonian for the fields on the circle corresponding to (Bf) as 



1 - r 2n ( 1 - \ 

H = rwf E N* + a E / dOxHe) -id 6 + - E KfiNp 

Q' a J ° \ Z (3 J 



xM, (12) 



with the representation of the particle numbers as 



N a = r d9xi{9)x«{0) ■ (13) 



o 

The equations of motion 

d tXa (9) = i [H, Xa {0)} = -a\d + i y £ \ a0 Nf)j Xa{0) , (14) 
on introducing new fields by 

M0)=e iE » Xaf, * f,0 Xa(9), (15) 



take a free form, 

(d t + ad e )M0) = 0. (16) 

To derive the commutation relations for the ip operators, it is convenient to use the 
identity 

e^^^x-y = X7 e^ Aa/3( ^~ W . (17) 
To prove this, we first note that the operators Xa and Xa satisfy 

[N a , Xp) = -$apXf3 , [N a , Xl) = -SapXp ■ ( 18 ) 

The first relation in (|18D written as N a Xp = Xp{N a — 6 a p), is generalized to N%xp = 
Xp(N a — S a /3) k , with k a positive integer, which straightforwardly leads to fll7|) . 



Using the identity (|T7D and similar identities for the operators xL we obtain from QTTj) 
the commutation relations for the ip fields: 

+ e- iX ^ d -^4(9')£(9) = (19) 
M0)4( 6 ') + e a ^ e - 6,) 4(6')M0) = S a pA a (9 - 9') , 
where the operators 

A Q (6 - 6') = e^e Xa ^ {e - e,) 5+ cv (6 - 6') (20) 
behave to some extent like (^-functions on the circle: 

2W d9'^ a {9')A a {9' -6) = tfjfi) , r de '^ e ~ = Md) • (21) 



/o Jo 
Consider the "Fourier transform" of the ip operators: 

^W = 4=EC^- (22) 

V Z7T K 

The (nonnegative) numbers n a are not integers, in contrast to the numbers n a in the 
expansion (fT0|). Allowed values of k q can be obtained by calculating the matrix elements 



of the field ip a between states with fixed particle numbers {n'p}. The only nonzero matrix 
elements correspond to the transition n' a + 1 — > n' a , without changing the other numbers 
of particles. Comparing the matrix elements of ip a , in its two representations, ( p2|) and 
QT5D with flrC|), yields the restrictions 

n a = n a + J2 Kpn'p , (23) 

/3 

where the nonnegative integers n a are the same as in (|10|). 

Let us introduce a set of basis states and label for a moment all the particles by a 
single sequence 1, . . . , N 

K N N ,---^T) = ^< N N ] ---< 1 M- (24) 

Then we obtain the restriction 

< = < + J2^Nf, i = l,2,...,N a , (25) 

where for each species a, % numbers particles in order of their appearance in fl24]), n l a 
correspond to the fermionic angular momentum quantum numbers (see (0)). We have 
also denoted by the number of creation operators for species (3 in the sequence fl2"l| ) 
(counted from the right to the left) standing before the operator creating the i th particle 
of species a (in other words, iVi is the number of particles of species (3 created before 
the i th particle of species a). 

Inserting the expansion ( ^2|) into (|19D, we obtain the commutation relations for the 
operators a" and (here we drop the index of k labelling the species of particles) 

aX + al Kp a: +Xap = 0, 

where II" are the projections on the subspace with the particle numbers {n'p} determined 
by (^3|). The projection operators satisfy the relations 
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For rational \ a p with common denominator q, explicit expressions for these operators can 
be derived: 



It is easily seen from (p8|) that on the states with fixed particle numbers, the operators 
(EH) have eigenvalues 1 if the relations (E^) are fulfilled and eigenvalues otherwise as it 
should be. 

The occupation number description determined by the creation and annihilation oper- 
ators and a" is not unique: as can be seen from (^), the interchange of any two of the 
creation operators leads to an equivalent state corresponding to the same fermionic quan- 
tum numbers {n l a }. Rearranging an iV-particle state does not change the total energy of 
the state, which is given by the sum over single-particle energies, 

£ = EE(K cff + ^)- (29) 

a j=l 

An ordering procedure may be introduced by demanding that each particle added 
into the system has to have a higher (or equal) value of k than the largest k already 
present in the system. This means that only the states (^4j) corresponding to the ordering 
K\ < ■ ■ ■ < are regarded as physical states, which makes the occupation number 
picture unique. Any of such normal ordered states can also be constructed in a unique 
way starting from a given set of the fermionic angular momentum quantum numbers {n l a } 
by putting particles into the system one by one and always making sure to put in the 
particle which gets the lowest possible value of k. 

Eq. fl25p then admits the interpretation that the spectrum of available single particle 
energy levels for a particle of a given species, say a, is shifted upwards by the presence of 
all particles below it. Each particle, say of species /3, causes a shift \ a p. By convention, 
particles of distinct species in the same level shift each other by \\ a p{— \\pa)- This means 
that the order of filling of levels of the same energy with particles of distinct species is 
unessential. 

Figure 1 illustrates the above normal ordering procedure in the case of a simple two- 
component system with integer statistics parameters. 
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3. Relation to exclusion statistics 

3.1. Thermodynamic limit of the free particle formulation 

With the ordering procedure described at the end of the previous section, Eq. fl25|) 
can be written in a compact form as 

^ = ^ + EE A ^K-^f), (so) 

where h(x) is the step function, which is defined here as h(x) = 1 for x > 0, h(x) = ~ for 
x = 0, and h(x) = for x < 0. We assume the ordering < < ■ ■ ■ < within each 
species (the resulting numbers nf within the same species are distinct). In the form (|3~0"D, 
the numbers nf are similar to "renormalized" quantum numbers introduced for integrable 



models of a Calogero type [|19 



We now specify the length of the circle as L and use the identification QTJ] 

a = ^v. (31) 

If one assumes that the harmonic potential is created in an anyon droplet of radius L/2n 
by an electric field, then the velocity v can be interpreted as the drift velocity on the edge 
E/B (the velocity of the edge excitations) where E = (m/e)u 2 R is the electric field on 
the edge. The thermodynamic limit is understood as uo — > and L — ► oo while v is kept 
fixed. We also introduce the momenta and pseudomomenta 

P-=fS Q , k * = T^ (32) 

distributed in the thermodynamic limit with the densities v a (p a ) and p a (k a ), respec- 
tively, so that the number of particles of species a in the interval (p a ,p a + dp a ) or in the 
corresponding interval (k a ,k a + dk a ) isQ 

v a (p a )dp a = p a (k a )dk a (33) 



1 The transition to the thermodynamic limit below is similar to that given in Ref. |2(|. We stress 
however that our starting point is the free field occupation number picture rather than the spectrum of 
an integrable model as in the above reference. 
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In terms of ©, Eq. (|30|) reads 



L & i 

or, in the thermodynamic limit, 

/■oo 

p a = k° - V / A Q/3 /l(fc a - k f3 ) P g(k ld )dk /3 . 

From this equation, we get 

dp a /dk a = l-J2K P Pp(k a ) 



(34) 



(35) 



(36) 



and with the use of (|3 



(37) 



1 — ^a/3Pf3 

The fermionic description of the energy levels implies that the (non-equilibrium) entropy 
has the form 

L 



S = — V / [-v a In v a - (1 - u a ) ln(l - u a )]dp a 
2ir „ Jo 



(38) 



or, in terms of p 



L ■ Z" 00 r r I r 

S = TtY, \-p a 1np a + l-J2pKpPi3 In l-E/3Aa/3P/3 

- [l - E/3(<W + Kp)pp] hi [l - E/3(<W + A a/3 )p /3 ] } d/c Q 

The latter expression, along with the expression for the energy 

L r°° 

E = 7rT, {\u c + vk a )p a {k a )dk a . 
2ir - JO 



(39) 



(40) 



shows that in the thermodynamic limit the model on the circle is equivalent to a gas 
of free chiral particles with linear dispersion (with a gap), the same for all the species, 
obeying ideal fractional exclusion statistics P, [5] with the fermionic statistics parameters 
We refer to the bosonic and fermionic statistics parameters |2l| g a p and \ a p as those 



appearing in the bosonic and fermionic counting of the statistical weight 



w = u 



(D b a + N a -1)\ 



n 



D f \ 



NJ(Db a -l)\ ^N a \(D f a -N a )\' 



(41) 
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with 



D b a = G a - J2 9^, Di = G a -J2 KpNp, (42) 
P P 

resulting in the expression for the entropy upon taking the thermodynamic limit 

N a — > oo, G a — > oo, with p a = N a /G a kept constant. Eq. (|4lD implies the relation 

between the bosonic and fermionic statistics parameters 

9a0 = Sap + . (43) 

Equations for the distribution functions p a in the equilibrium can be derived from 
Eqs. (p9|,^D in the usual way. In terms of the single state grand partition functions £ Q , 
related to the p a 's by || 

ta = l + z v ,/° Vn , (44) 

these read 

(&* - 1) II & 8 = e^-^- vka)tT . (45) 
P 

3.2. Heat capacity 

We now assume that all the species are kept at the same chemical potential, \x a = p. 
We also restrict to the case where all the distributions p a (and, correspondingly, £ Q ) are 
equal to each other. It then follows from (|45| ) that 

M 

\M = Y, Kp (46) 
p=i 

is independent of a, and ([45]) coincides with the equation for the partition function for a 
single species of particles with the statistics parameter (pEE|): 

(6a - l)6a M = e { ^-^- vk)/T . (47) 
We calculate the low temperature heat capacity for this case following essentially 



Ref. p2[ for the case of a constant density of states, which holds for a linear dispersion in 
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one dimension (the difference from Ref. |22j is that in the case at hand the particles are 
chiral and have a gap in the dispersion law). Introducing the variable w by 

1 



1 — w 



(48) 



we represent the particle density d = N/L = J2 a Jo° Pa(k)j^ and the energy (fjQ|) as 
integrals over the variable w: 

MT r«(o) dw 



d 

2irv Jo 1 — w 
E MT r°(P) dw 



L 2ttv Jo 1 — w 



\p a + (1 + AM)Tln(l -w)- Tlnw] 



(49) 



where w(0) is the solution of Eqs. (47) and (^q) for k — 0. From the first equation in (^9| 
we find 



w(0) = 1 - e 



-2-wvd/MT 



(50) 



From Eqs. ([47D , (|48 ) we then have 



/x a = icj c + 2vr(AM + l)v— + T In (l - e 



-2-wvd/MT 



Finally, combining this with the second equation in (fi9|), we obtain 



E 

L =d 



\u c + (AM + l)vrt;^: 



m „ , x MT 2 Mo) , lnw 
+ Tdlnw(0) - / dw- 



2rcv Jo 



1 — w 



(51) 



(52) 



At low temperatures T 2nvd/M, up to non-perturbative corrections, containing the 
exponential factor e - 27Tvd / MT ; ^(0) can be set to one. The integral in (52) then equals 



— ^7r 2 , which results in the heat capacity 



L Qv 



(53) 



In the single species case this result reduces to that obtained in Ref. in a different way. 
The present derivation shows that (|53D exhausts all the perturbative (behaving as powers 
of the temperature) terms in the heat capacity. 



3.3. Asymptotic Bethe ansatz equations 
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The connection with ideal exclusions statistics suggests that the equations of the model 
can be represented in a form similar to asymptotic Bethe ansatz equations |L6| (the latter 
have the form of the thermodynamic Bethe ansatz equations |23[, with two-body scattering 
phases (matrices) having a step- wise form). Indeed, introducing the "dressed" energies 
e a by 

1 



a e^l T + 1 

so that 



(54) 



£ a = 1 + e~ ealT , (55) 

we find that Eqs. (|4"5| ) become 



e a = -Ha + e° a (k a ) In (l + e^-^ )/T ) , (56) 

where 

e° a {k a ) = \uj c + vk a . (57) 
Equations (|56|) , along with (|35|) , have the form of the asymptotic Bethe ansatz equa- 



tions for a system of particles with the "bare" energy ( 57]) and the step-wise "two-body 
scattering phases" 

0(kj - fcf ) = Kph(kf - fcf ) . (58) 
One can in addition introduce the density of holes p^(/c), related to the interval dk a , by 

(l-u a ( V a ))d V a = p h a (k a )dk a . (59) 
With the use of (|33|) and (pB]), we obtain 

^i-EfV + M^. (so) 

For T = 0, in the case of equal boundary pseudomomenta, Uq = k , we have 

Pa = E/3(-W + Sap)' 1 , Pa = 0] k < k ] 

, (61) 
Pa = 0, p£ = 1; k> k . 
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We now turn again to the above special case of equal p a for all species and evaluate 
the Fermi velocities of quasiparticles (the same for all species) in the picture of interacting 
fermions for T = 



Vp = 



de a 



dp 



P=Pf 



dk de a 
dp dk 



(62) 

k=kg 



where pp is the Fermi momentum for the distributions v a and k§ is the boundary (pseudo 
Fermi) momentum for the distributions p a . 

For zero temperature, because of the jumps in the distributions of pseudoparticles and 
holes, the expression ( |5"2"D is not well defined for k = k$. It is well defined however for 
k = kg(= ko — 0) and for k = k^(= k + 0). In the above special case of equal distribution 
functions for all the species, we obtain from (|56"D, s' a {kQ ) = v(\M + l) -1 and s' a {kQ) = v. 
With dp/dk = (AM + l) -1 for k = k$ and dp/dk = v for k = k$ following from (^), we 
conclude that 

Vf{K) = v f(K) = ^ for T = 0. (63) 
Representing then the Fermi velocities fl62|) as 



k? = { 4W/pW. * = Af, (64) 

£ a (k), k = k , 

we remark that the first line in (pip, which is normally used to evaluate the Fermi velocities 
in the thermodynamic Bethe ansatz method (see e. g. 11241), in the case of step- wise two 
body phase shifts (asymptotic Bethe ansatz) holds only for k < k$. 

Similar peculiarities due to a step-wise form of the two-body scattering phase also arise 
for other quantities determined by the asymptotic Bethe ansatz equations. We discuss here 
the dressed charge matrix (which is related to Friedel oscillations, conductivity etc., see 
p5|). With the two-body scattering phases (|58D, the dressed charge matrix is determined 
by the equations [p4| 

rko 



Zapikp) = 5 a/3 - / Za 7 (& 7 ) \ 1 p8{k 1 - kp) dk r (65) 
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The function Z a p(kp) thus has a jump at the boundary pseudomomenta k Q : below the 
pseudo Fermi level 

Za(3(K ) = (fi ,_1 )"/3 5 (66) 

while above the pseudo Fermi level 

Z a p(k^) = 5 a/3 . (67) 

The relation between the dressed charge matrix and exclusion statistics matrix can also 
be written in a matrix form as 

Z{k»)Z{k+) = g- 1 . (68) 
which is different form the result Z(k )Z(k ) = g _1 obtained in Ref. . 



Note that the representation of the Fermi velocity ( p3|) and the relation of the dressed 
charge matrix to the statistics matrix ( |68|) are valid for an arbitrary "bare" energy of 
particles E° a {k) (which in the case at hand is fl57|)). For the case of a single species, the 
relation similar to (|68"D, between the dressed charge function and the fractional exclusion 
statistics parameter, was discussed in Ref. p7 |. 



3.4- Thermal excitations 

Because of the nontrivial temperature dependence of the dressed energy e(k) deter- 
mined by Eq. fl56|) , it is not obvious that the above zero temperature quasiparticle exci- 
tations survive at finite temperatures. In this subsection, using the exclusion statistics 



representation of the model, we calculate the Fermi velocities fl6"2|) at low but nonzero 
temperatures. The derivatives with respect to k in ( j62| ) then have to be evaluated at the 
point fi a = £° a {k)- 

The derivative dp/dk is given by (p6[), which in the case at hand reads 

j| = 1 - XMp a . (69) 

To calculate the derivative de a /dk, we note that, according to (P5|) and (|55|), e a depend 

on k only via x = e^ Q_e °^ fc ^ //T ', which yields 

de a _ dxcka _ yx ^a 1 ^ 
dk dk dx dx ^ a — 1 
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Taking into account the relation between the single state distribution functions and par- 
tition functions (£ = Yia^a), 

p = x-^ = Mx-^ (71) 

and expressing £ Q in terms of p a with the use of we get 

de 

" u(l-AMp a ). (72) 



dk 

Note that both the derivatives fl69|) and (|72|) are well defined at the point fi a = e° a {k). 
Using these, we finally obtain the Fermi velocities of quasiparticles for T^O 

v a F = v . (73) 

This calculation proves that the Fermi velocities of the quasiparticles at nonzero temper- 
atures remain the same as at T = 0. As the derivation shows, the conclusion holds for an 
arbitrary bare energy e Q a {k) in the asymptotic Bethe ansatz equations (|56|) . 

4. Bosonization and relation to edge excitations in multilayer FHQ states 

In this section we bosonize the model governed by the Hamiltonian (0) generalizing 
the procedure given in Ref. Jl| to the case of several chiral fields. We then discuss the 
relation of the bosonized form of the model to the effective low energy description of edge 



excitations in multilayer FQH states [17, 



Consider the Hamiltonian (|T^) of chiral fields on the circle of length L specified by 
(pTTI). The fact that the model is described with only one velocity v means that in the 



picture of the anyon droplet, all the species of particles have the same boundary radius. 
We assume that the boundary (pseudo Fermi) values of k in the ground state are the 
same for all species. It then follows from the occupation rules and the ordering procedure 
described in Section 2 (or, equivalently, from Eq. (|30|) ) that the pseudo Fermi angular 
momentum is 



/,;; - (X aa + l)(N 0a - 1) + (74) 



.a 
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where N 0a is the number of particles of species a in the ground state. We introduce the 
chemical potentials for zero temperature by 



with 



A*o = + aK F 



K a F = (X aa + l)[N 0a --)+ J2 KpN w . 



(75) 



(76) 



This choice corresponds to the chemical potential lying on the midway between the highest 
occupied and the lowest unoccupied energy levels. 

The pseudo Fermi levels in the model on the circle are directly related to the physical 
size of the anyon droplet: for large N a one has 



(77) 



where r = \j2/eB e s is the effective magnetic length. This means that all k f have the 
same value, Kp, from which follows that J2p 9ap^op should not depend on a. This also 
implies that //q = £*o- 

We fix the ground state subtracting /i N from the Hamiltonian ([121), 



-id e - k f ) + ]- 

1 p 



Xa(0) . (78) 



By normal ordering the Hamiltonian with respect to the Fermi levels, subtracting a con- 



stant and redefining the field operators Xa for the rest of this section by Xa — ► e 
we obtain the Hamiltonian describing the excitations 



-iN Oa 



^ d9 : xl(0) (-id e + ~) Xa(d) : +\ £ ^pQ»Qp 



Xai 



(79) 



where Q a = N a — N 0a are the charge operators. 

The fermion fields Xa admit the representation (for details of the bosonization proce- 
dure, see Ref. pSfl ) 



f//.\ e2 



xUo) 



(80) 
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where the chiral boson fields are given by 



n>0 



Cl an € 



in6 



u an e 



51) 



with [0o = id a@ an d Km,4 n ] = fiapfimn- The normal ordering in (|80|) refers only to 
a an and a^. The fields Q satisfy the commutation relations 



l«(0),0/j(0')] =^sgn per (^-^ 



(82) 



where sgn (#) is the periodic sign function: (d/dQ) sgn per (6*) = 2S per (9). The extra phase 
in flS0P containing the charge operators ensures that the fermion operators for distinct 
species anticommute (rather than commute) with each other. 
The bosonized form of the Hamiltonian (|79|) reads 



H' 



-E 



i 

47T 7o 



2- 



d0 : (dg 



an"'*'! 



13) 
»4) 



where p Qi g are the bosonic statistics parameters (see fl43f)). 

Since g Q /3 is a symmetric matrix, it can be diagonalized by an orthogonal transfor- 
mation, (0~ 1 gO) a p = A a 5 a p. The Hamiltonian then becomes that of a set of uncoupled 
components with (bosonic) statistics parameters A a , 



J5) 



where a an = 0~pap n , Q a = 0~pQp. The transformed operators a am and oL still obey 
the bosonic commutation relations. 

At this stage one has to demand that all the eigenvalues A a of the matrix g a p should 
be positive, in order for the Hamiltonian fl35|) to be positive definite. This means, as we 
will see below, that all the edge modes in our model propagate in the same direction. 
Note that the diagonal form of the Hamiltonian (|85|) implies that the low-temperature 
heat capacity is just a sum of heat capacities from all the species, nT/Qv from each, in 
agreement with (|53|) . 
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The Hamiltonian (|85| ) can be transformed to a free form 



u 2 



4tt 

where the new Bose fields are defined by 

$ a {e) = -^=4>% - y/T a Q a e + ^e4= \^ n6 ~ ~ a i^' in6 \ , (87) 

with [(f>%, Qp] = i5 aP . 

One can define the charged operators corresponding to the fields Q : 

These operators have charges q a = 7«/v / A« (determined by the commutator X«]) as 
well as "statistics" 7^ in the sense that 

xtmW) = e l ^ s ^ {d - e,) xU0lxU0) ■ (89) 

In particular, if one chooses 7^ = \/A^ an d if A Q = 2m a + 1, with m a a positive integer, 
the operator (RSI) is an "electron" operator of the kind discussed by Wen H. Another 



special case r y a = 1/ \fA^ =1/ \/2m a + 1 corresponds to a particle with fractional charge 
and statistics equal to l/(2m a + 1) like a fundamental quasiparticle in a v = l/{2m a + 1) 
FQH layer. The fields corresponding to distinct species anticommute with each other, 

A general charged operator can be composed of operators of the above type: 
which corresponds to the total charge J2 a p Oaplp/ associated with the operator 

Qa- 



To discuss the relation to the edge excitations in multilayer FQH systems $J\, §J, 
we rescale the Bose fields 4> a and transform back to the non-diagonal picture, using the 
inverse of the above orthogonal transformation. This defines new fields 

4> a = O af3 ^=~$p , (91) 
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obeying the commutation relations (cf. (B3)) 



iMM?)] =iirg- 1 pSgn pcr (9-9') 



(92) 



The Hamiltonian 



takes the form 



(93) 



For the case where the diagonal elements of the matrix g a p are odd numbers and 
the off-diagonal elements are integers, Eqs. (|92|), (|9~3"D reduce to those describing edge 
excitations in multilayer abelian FQH states H if one identifies 



9aP = K a/3 , V9a/3 = K 



a/3 , 



(94) 



where K a p is the topological matrix (in the symmetric basis) and V a p is a positive definite 
matrix describing the interaction between the chiral bosonic modes 0. Thus, our model 
with only one velocity v, corresponds to the special case where the matrices K and 
V are proportional to each other. In the more general case of the effective Luttinger 
liquid description of the edge excitations, when different edge modes may have different 
velocities, the matrices V and K can be simultaneously diagonalized, and the velocities 



are determined by the ratios of their respective eigenvalues [29 



Note that the conclusion about the identification of the exclusion statistics matrix 
with the topological matrix (the first relation in (0)) was previously obtained by Fukui 
and Kawakami for hierarchical FQH states (by comparing the excitation spectrum for the 
edge excitations and chiral particles obeying ideal fractional exclusion statistics) [T2"|] . 

In addition, expressing the general charged operator (|90|) in terms of the fields (fi a as 



x\e) 



(2tt) a/ /2 



(95) 



with l a = Qj a-y / /Ag7 / g, we find that this operator has the total charge 



Ql = l 09al 



a/3 



(96) 
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and statistics Y^ a p^ag a 8h- For integer l a , this agrees with the expressions obtained by 
Wen and Zee in the symmetric basis fL%. Finally, the filling factor is 



f = E^ 1 = E^T L *« ( 97 ) 

a/3 a A a 

where t a = Y^pO~\. Figure 1 illustrates the first equality in fl9"7D for a simple example of 
two components. 

5. Concluding Remarks 

We have constructed a model describing several chiral fields on the circle which has 
the same quantum numbers as a system of anyons of several species in the lowest Landau 
level. The parameters of the model range in wider intervals than the original (statis- 
tics) parameters of anyons; this is achieved by the analytic continuation of the solutions 
corresponding to anyons in the LLL. 

The model incorporates some features of the physics described by the electron wave 
functions corresponding to abelian multilayer FQH states. The harmonic potential added 
into the system of anyons plays the role of the confining potential of electrons in the FQH 
states. The model on the circle was found to recover the effective low energy description of 
edge excitations in the special case of equal velocities of all edge modes. In this sense, the 
model obtained can be considered as a possible dynamic theory underlying the effective 
chiral Luttinger liquid description of edge excitations. 

In this context we note that some of the correlation functions as well as low tem- 
perature thermodynamics of edge excitations have been obtained using only the effective 
Luttinger liquid description of the edge excitations (see e. g. the review |4j]). On the other 
hand, recent calculations of nonequilibrium transport properties through a point contact 
in a Luttinger liquid (which is related to tunneling transport between edge states in FQH 
devices) are based on a particular (integrable) dynamic model of the Luttinger liquid [|30| 



In this connection the question arises to what extent these nonequilibrium properties 
depend on a particular dynamic model of edge excitations. The model we have discussed 
has a simple dynamics encoded in step-wise two-body scattering phases when it is for- 
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mulated in the form of the asymptotic Bethe ansatz equations (see Sect. 3). It seems 
therefore to be interesting to use this model to investigate the above issue. 

Another remark concerns the edge excitations for hierarchical FQH states. Integrable 
models with long range interactions can be constructed which are described by the same 
matrices as the topological matrices corresponding to hierarchical FQH states (see ||31|| ). 
A simple link between the integrable models with long range interactions and edge ex- 
citations for FQH states is given by the (ideal) fractional exclusion statistics. Exclusion 
statistics provides, in addition, a connection with anyons in the LLL ||, |l^, |5| . This sug- 
gests that an explicit dynamic model of edge excitations for the hierarchical FQH states 
can also be constructed, along the lines of the present paper, starting from the appropriate 
picture of anyons in the LLL. 
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Fig. 1: Ground state of a two-component system with gu = 3, #22 = 5 and gu = 
#21 = 1 (bosonic representation). The white and black circles denote component 1 and 
component 2, respectively. The state has been built up according to the normal ordering 
convention, with increasing energy (k) values. The pattern repeats periodically. The 
filling factor is 3/7 = J2 a p9a8 the number of particles in the ground state is such that 
the Fermi level is the same for both species). 
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